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Abstract 

Experimental data are reported on moisture diffusion and the elastoplastic re- 
sponse of an intercalated nanocomposite with vinyl ester resin matrix and mont- 
morillonite clay filler at room temperature. Observations in diffusion tests show 
that water transport in the neat resin is Fickian, whereas it becomes anomalous 
(non-Fickian) with the growth of the clay content. This transition is attributed 
to immobilization of penetrant molecules on the surfaces of hydrophilic clay layers. 
Observations in uniaxial tensile tests demonstrate that the response of vinyl ester 
resin is strongly elastoplastic, whereas an increase in the clay content results in a 
severe decrease of plastic strains observed as a noticeable reduction of curvatures of 
the stress-strain diagrams. This is explained by slowing down of molecular mobility 
in the host matrix driven by confinement of chains in galleries between platelets. 
Constitutive equations are developed for the anomalous moisture diffusion through 
and the elastoplastic behavior of a nanocomposite. Adjustable parameters in these 
relations are found by fitting the experimental data. Fair agreement is demonstrated 
between the observations and the results of numerical simulation. A striking simi- 
larity is revealed between changes in diffusivity, ultimate water uptake and the rate 
of plastic flow with an increase in the clay content. 
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1 Introduction 



This paper is concerned with modelling moisture diffusion through a hybrid compos- 
ite consisting of a polymer matrix filled with clay nanoparticles. Permeation of small 
molecules through polymer-clay nanocomposites has attracted substantial attention in 
the past decade, see, e.g., [I], ||, |3|, f|, |5|, ||, [7], ^ |H], [ll], [121 , to mention a few. This 



may be explained by a strongly increasing interest to development and characterization of 
nanocomposites whose physical properties demonstrate a noticeably improvement com- 
pared to those of pristine polymers at rather low (conventionally, less than 5%) filler 
concentrations. 

In the previous works, barrier properties of polymer-clay nanocomposites were mea- 
sured by the permeability cup method |1|, f|, ^ [T(J and the microgravimetric technique 
|| [Ll], |I2| . The former approach characterizes the steady-state mass transport only. In 
the latter case, the microgravimetric analysis focused only on (i) the initial parts of the 
mass gain curves (which provided estimates of diffusivity for the Fickian diffusion) and 
(ii) the equilibrium water uptakes. This implies that the kinetics of moisture diffusion in 
polymer-clay nanocomposites have not yet been investigated in detail. 

Changes in diffusivity and permeability of nanocomposites with clay content are con- 
ventionally explained within the concept of tortuous paths According to this ap- 



proach, the path that a small molecule of a penetrant must travel in a polymeric matrix 
substantially increases in the presence of intercalated and/or exfoliated clay layers (that 
possess aspect ratios of order of 10 3 ). This results in a noticeable decrease in the coefficient 
of diffusion (estimated as an average value of the square of the end-to-end distance passed 
per unit time). According to the Nielsen equation |13| , the tortuosity factor (the ratio of 
diffusivity of a composite to that of the pristine polymer) is inversely proportional to a 
linear function of the volume fraction of filler. Among other phenomenological relations, 
it is worth noting the Cussler formula |14| (the tortuosity factor is inversely proportional 
to a linear function of the square of the volume fraction of filler), the Barrer formula 
(the tortuosity factor is proportional to the square of the volume fraction of polymer), and 
the power-law equation (the tortuosity factor is proportional to a fractional power of 
the volume fraction of matrix). The Nielsen formula was recently used to fit experimental 
data for polyimide-clay || and polyester-clay || nanocomposites. The Barrer equation 
was employed to match observations for Nylon 11 filled with nanosize silica ||. In all 
cases, the quality of approximation was far from being excellent. 

To improve the accuracy of fitting, Fredrickson and Bicerano [pl| and Bharadwaj 
17j proposed phenomenological relations that account for the orientation of platelets, 



their shape, and the state of delamination. This (merely geometrical) approach may 
reduce deviations between experimental data and results of numerical simulation (due 
to an increase in the number of adjustable parameters), but it cannot overcome the 
main shortcoming of the original concept: the tortuous path theory is grounded on the 
assumption that the presence of nanoparticles does not affect the diffusivity of a polymer 



matrix 16 



This postulate appears to be rather questionable, because observations show that 
mass transport in polymers is intimately connected with their viscoelastic properties |18|] . 
Experimental data in oscillation tests with small strains reveal an increase in the storage 
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modulus and a decrease in the loss modulus with filler content, v, for nanocomposites 
with Nylon 6 ||19|| , Nylon 11 ||, Nylon 1012 [l^| and polyester || matrices. The reduction 



of the loss tangent (that reflects changes in the viscoelastic response) is accompanied 
by a decay in viscoplastic properties revealed as formation of a pronounced yield point 
for polypropylene-clay nanocomposites and transition from ductile to brittle fraction in 



Nylon 6-clay nanocomposites with the growth of the filler fraction [113. These observations 
demonstrate that mobility of macromolecules in a polymeric matrix slows down with the 
clay content. Adopting the concept of "chain immobilization factor," see |DJ and the 



references therein, one can conclude that this deceleration of molecular mobility should 
be accompanied by a decrease in diffusivity of small molecules, which is not accounted 
for by the concept of tortuous paths. 

The objective of this study is three-fold: 

1. To report experimental data in moisture diffusion tests and in uniaxial tensile tests 
on vinyl ester resin-montmorillonite clay nanocomposite with various contents of 
filler. 

2. To derive constitutive equations for water sorption and stress-strain relations for 
uniaxial deformation of a nanocomposite and to find adjustable parameters in the 
governing equations by fitting the observations. 

3. To demonstrate that the clay fraction affects the coefficient of diffusion and the 
constants responsible for the viscoplastic behavior in a similar way. 

The choice of vinyl ester resin for the experimental analysis is explained by numer- 
ous applications of this thermosetting polymer as a matrix for glass-reinforced polymer 



composites employed in construction and repair of civil structures |20 . 

Montmorillonite (MMT) is an inorganic clay conventionally used for preparation of 
hybrid nanocomposites. It possesses a layered structure constructed of two tetrahedral 
sheets of silica surrounding an octahedral sheet of alumina or magnesia. The layers (with 
thickness of 1 nm) are stacked by weak dipole forces, while the galleries between the layers 
are occupied by metal cations. 

To evaluate the effect of clay content on the mobility of macromolecules in the matrix, 
we concentrate of water sorption tests and uniaxial tensile tests. Our experimental data 
(Figures 1 to 5) demonstrate that moisture diffusion in the neat vinyl ester resin is Fickian. 
With an increase in the clay concentration, the Fickian diffusion is transformed into 
an anomalous transport of the penetrant molecules. Anomalous diffusion of water in 
polymers and polymer composites has been recently studied in [^1], [23|, |24|, |25| . 



We suppose that this anomalous moisture uptake in a nanocomposite filled with MMT 
particles may be explained by immobilization of water molecules on the surfaces of hy- 
drophilic clay layers. A model for a two-stage diffusion with immobilization of penetrant 
molecules was proposed by Carter and Kibler [^6j and Gurtin and Yatomi ||27|| . Unlike 



these works, we assume that bounded molecules cannot leave the sites where they were 
immobilized, which allows the number of adjustable parameters to be substantially re- 
duced. The diffusion process is described by differential equations with three material 
constants that are found by matching the water uptake curves plotted in Figures 1 to 5. 
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Our observations in tensile tests (Figures 6 to 10) reveal that the stress-strain curve 
for a pristine polymer substantially differs from a straight line, which means that the 
mechanical response of the neat vinyl ester resin is strongly elastoplastic. With an increase 
in the clay content, the stress-strain diagrams become "less curved," which is associated 
with a decay in the material plasticity. Similar observations were recently reported for 
nanocomposites with Nylon 6, polypropylene [I9~|, and poly(e-caprolactone) || matrices 
filled with MMT clay. 

To describe the elastoplastic behavior, we treat a nanocomposite as an equivalent 
network of chains bridged by junctions (entanglements, physical and chemical cross-links, 
and clay platelets). The elastic response is attributed to elongations of strands in the 
network, whereas the plastic response is attributed to sliding of junctions with respect 
to their reference positions in the medium (the non-affine deformation of the network). 
The decrease in the curvatures of the stress-strain diagrams with the clay content is 
ascribed to a decay in mobility of junctions whose motion is restricted by the presence of 
nanoparticles. A model is developed to describe constrains imposed by the clay platelets 
on the mobility of chains in the host matrix. The constitutive equations are determined 
by three adjustable parameters with a transparent physical meaning. These parameters 
are found by matching the stress-strain curves in uniaxial tensile tests with small strains. 

Unlike the previous study |28| , where observations were reported for a nanocomposite 
filled with MMT clay treated with dimethyl benzyl hydrogenated tallow quaternary am- 
monium chloride (which implied that the clay layers were not chemicaly bonded to chains 
in the host matrix), the present work focuses on transport of water molecules through a 
nanocomposite filled with MMT clay modified with vinyl benzyl ammonium chloride (the 
latter ensures that some macromolecules are chemically bonded to the clay platelets). 

The exposition is organized as follows. Observations in moisture diffusion tests and 
uniaxial tensile tests are reported in Section 2. A model for water diffusion through a 
nanocomposite is derived in Section 3 and its material constants are found in Section 4. 
Constitutive equations for the elastoplastic response of a nanocomposite are developed in 
Section 5. Adjustable parameters in the stress-strain relations are determined in Section 
6. A brief discussion of our findings is presented in Section 7. Some concluding remarks 
are formulated in Section 8. 



2 Experimental procedure 
2.1 Preparation of samples 

The polymer used was DERAKANE 411-350 epoxy vinyl ester resin (Dow Chemical 
Co.) containing 45 wt.-% of dissolved styrene. To cure the resin at room temperature, 
it was mixed with 0.5 wt.-% of 6% cobalt naphthenate catalyst (Sigma Aldrich Co.). 
Additionally, 0.05% of 99% N,N dimethyl aniline (Lancaster Synthesis, Pelham, NH, 
USA) was used as an accelerator, and 1.5% of methyl ethyl ketone peroxide with 9% of 
active oxygen (Sigma Aldrich Co.) as an initiator. 

A sodium montmorillonite clay Cloisite NA (Southern Clay Products Inc., Gonzales, 
TX, USA) was treated with vinyl benzyl ammonium chloride (Sigma Aldrich Co.). 30 g of 
Cloisite NA was dispersed in 800 ml of distilled water. Separately, 6.35 g of vinyl benzyl 
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ammonium chloride (that corresponds to 100 meq/100 g of clay) was dissolved in 80 ml 
of distilled water. After stirring of the clay dispersion with a mechanical stirrer for 5 h, 
the vinyl benzyl ammonium chloride solution was added drop by drop in the dispersion. 
The mixture was stirred for another 3 h, and the organo-MMT clay was filtered using a 
vacuum filter, dried in an oven at 100 °C, and ground in a ball mill. 

Neat resin coupons were cast by pouring the reaction mixture into Teflon molds with 
the dimensions 50 x 12.5 mm and the thickness ranging from 0.2 to 0.6 mm. The organ- 
ically treated clay was added to the liquid resin and manually stirred. The mixture was 
then degassed in a vacuum oven to remove air bubbles. Afterwards, the catalyst, initiator 
and accelerator were added. The mixture was allowed to cure at room temperature for 
24 hours, and it was subsequently post-cured in an oven for 3 hours at 90 °C. 

TEM (transmission electron microscopy) micrographs of samples show that an inter- 
calated nanocomposite is produced |29[. It is worth noting, however, that the thickness 
of interlayers between platelets does not increase substantially, which may be explained 
by cross-linking between unsaturated sites on the surfaces of clay layers and vinyl groups 



of chains in the host matrix [29, 30 



2.2 Diffusion tests 

Diffusion tests were performed by immersing samples with rectangular cross-section in 
distilled water at room temperature (T = 25 °C). The samples had a dry mass ranging 
from 120 to 400 mg. They were stored in a controlled humidity chamber, and contained 
0.05 ± 0.005 wt.-% of water at the beginning of diffusion tests. The samples were period- 
ically removed, blotted dry with a lint-free tissue, weighed and re-immersed into water. A 
typical experiment lasted ten days, and, on the first day, readings were taken as frequently 
as every 15 min. The balance used had an accuracy of 1 mg, and 3 to 5 replicate runs 
were carried out for a given set of condition. 

Experimental data are presented in the form of a relative mass gain, $, as a function 
of a reduced time, t. The reduced time, i, is given by 

t , s 

where t is time elapsed from the beginning of a test, and 21 is thickness of a sample. The 
relative moisture uptake, reads 

K[oo) 

where R(t) is the mass gain at the reduced time t and R(oo) is the maximal mass gain. 

Typical moisture uptake curves are plotted in Figures 1 to 5. These figures demon- 
strate that diffusion of water molecules in the neat resin is Fickian, but it becomes strongly 
anomalous with an increase in the clay content. 

2.3 Mechanical tests 

Dumbbell specimens for mechanical tests were prepared in accordance with ASTM D638 
specification. Tensile tests were performed at ambient temperature by using a 100 kN 
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Instron machine model 8501 with a cross-head speed of 0.254 mm/min. This strain rate 
ensures practically isothermal loading conditions. The strain was measured independently 
using a strain gauge affixed to the mid-point of a specimen. The tensile force was measured 
by a standard load cell. The engineering stress, a, was determined as the ratio of the 
tensile force to the cross-sectional area of a specimen in the stress-free state. For any 
clay concentration, u, four specimens were tested. Typical stress-strain diagrams are 
presented in Figures 6 to 10. These figures show that the mechanical response of the neat 
resin is strongly nonlinear, but the nonlinearity of the stress-strain curves is reduced with 
an increase in the clay content. 

DSC (differential scanning calorimetry) measurements show that the glass transition 
temperature of the neat vinyl ester resin is about 98 °C. It increases with the clay content 
reaching 112 °C at v = 1.0 wt.-%, and remains practically constant at higher concen- 
trations of filler P5| . At room temperature, the nanocomposite is far below its glass 
transition point, which implies that we can disregard its viscoelastic response and at- 
tribute the nonlinearity of the stress-strain curves to the elastoplastic behavior. 



3 A model for moisture diffusion 

A sample is treated as a rectilinear plate with thickness 21. We introduce Cartesian 
coordinates {x, y, z}, where the axis x is perpendicular to the middle plane of the plate, 
and the axes y and z are located in the middle plane. Length and width of the plate 
substantially exceed its thickness, which implies that the moisture concentration depends 
on the only coordinate x. 

After the plate is immersed into water, three processes occur in the nanocomposite: 

1. sorption of water molecules to the sample faces from the surrounding, 

2. diffusion of the penetrant into the plate, 

3. adsorption of water molecules on the hydrophilic surfaces of clay layers, where these 
molecules become immobilized. 

It is conventionally accepted that the rate of sorption in glassy polymers noticeably 
exceeds the rate of diffusion |[24|| , which implies that the sorption equilibrium is rapidly 
established. The equilibrium condition reads 

n(t,x) x=± = n°, (3) 

where n is the moisture concentration at time t at point x [the number of water molecules 
in the polymeric matrix occupying a volume with the unit area in the plane (y, z) and 
the thickness dx reads n(t, x)dx ], and n° is the equilibrium moisture concentration in the 
matrix on the faces of a sample. 

Diffusion of a penetrant through a matrix is described by the mass conservation law 

dn dJ drii 

~dt = ~dx~~ In' ^ 
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where rii(t,x) denotes the concentration of water molecules immobilized at the surfaces 
of clay layers, and the mass flux, J(t,x), is given by the Fick formula 

J = D w < 5 > 

The constant D in Eq. (5) stands for the coefficient of diffusion through the polymer 
(which means that only the mass transport through the host matrix is taken into account). 
It is assumed, however, that D depends on the clay content: the presence of clay platelets 
in the polymeric matrix leads to a decrease in the chain mobility that provides a driving 
force for diffusion. 

Adsorption of water molecules on the hydrophilic surfaces of clay layers is governed 
by the first-order differential equation 

dni 



dt 



kn{n\ — ni). (6) 



Equation (6) means that the rate of immobilization of small molecules is proportional 
to (i) the penetrant concentration in the matrix, n, and (ii) the number of "unoccupied 
sites" on the surfaces of clay platelets, n\ — ni, where n\ stands for the total number of 
sites where water molecules are bounded. The rate of adsorption, k, and the maximal 
concentration of unoccupied sites, n°, are assumed to depent on the clay content, because 
the presence of nanoparticles affects mobility of chains in the polymeric matrix, and, as 
a consequence, its chemical potential. 

The following initial conditions are accepted: 

n(t,x) =0, ni(t,x) =0, (7) 



t=o 



t=o 



which mean that the moisture content in a sample before testing is disregarded. 
The water mass gain (per unit mass of the sample) is given by 



Ki ['^ r 

R(t) = — n(t,x) +n 1 (t,x) dx, (8) 

Zil J — / 



21 p J-i 

where k is the average mass of a penetrant molecule, p is mass density of the polymeric 
matrix, and the functions n and n\ obey Eqs. (3) to (7). 

Our aim now is to transform the governing equations. Introducing the dimensionless 
variables 

x n rt\ 

X = — , C = , Ci — — 

I n° n° 

and combining Eqs. (1), (4) and (5), we arrive at the mass flux equation 

dc _ (?C _ dCi 

dt dx* at' [) 



where 

In the new notation, Eq. (6) reads 



dci 
~dt 



D = AD. (10) 



Kc(C- Cl ), (11) 
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where 

K = 4kl 2 n°, C=X (12) 

It follows from Eq. (7) that the initial conditions for Eqs. (9) and (11) are given by 

c(0,x) = 0, ci(0,x) = 0. (13) 

Bearing in mind that the functions c(i,x) and ci(i,x) are even functions of x, we find 
from Eq. (3) that 

dc - - 

— (f,0)=0, c(M) = l. (14) 

Equation (8) implies that 



R(t) = -n° [ 1 \c(t,x) + c 1 (t,x)]dx. (15) 
p Jo J 



P 

According to Eqs. (9) and (11), with an increase in time, t, the functions c and C\ approach 
their steady-state values 

c(oo,x) — l, ci(oo,x) — C. (16) 
Equations (15) and (16) result in the formula for the maximal moisture uptake, 

R(oo) = -n°(l + C), (17) 
P 

The relative mass gain is found from Eqs. (2), (15) and (17), 

$ (^) = TT~7^ f 1 \c{t,x) + c 1 (i,x)]dx. (18) 
1 + G jo L J 

Equations (9), (11) and (18) with initial conditions (13) and boundary conditions (14) 
are determined by 3 adjustable parameters: 

1. the reduced diffusivity D, 

2. the rate of moisture adsorption, K, on the hydrophilic surfaces of clay layers, 

3. the reduced concentration of unoccupied sites, C, on the surfaces of clay platelets. 
These quantities are found by fitting the experimental data in diffusion tests. 

4 Fitting observations in diffusion tests 

We begin with matching the experimental data for the neat vinyl resin depicted in Figure 
1. In the absence of clay particles, one can set C = 0. This equality together with Eqs. 
(11) and (13) implies that Ci(t,x) = 0. Bearing in mind this formula, we find from Eq. 
(9) that the function c obeys the conventional diffusion equation 

dc ^<9 2 c . . 

at = D a¥ < 19 > 
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with initial condition (13) and boundary conditions (14). The initial-boundary problem 

(13) , (14) and (19) is solved numerically by the finite-difference method with an explicit 
algorithm. We divide the interval [0, 1] into M subintervals by points x m = mAx (m = 
0, 1, . . . , M) with Ax = 1/M, introduce discrete time t n = nAt, and replace Eqs. (13), 

(14) and (19) by their finite-difference approximation 



C\tn+li%m) c(t n , X m ) + D A ^ c(t n ,X m +i) 2c(t n ,X m ) 



c(0,x m ) = (m = 0,1,..., M-1), c(0,x M ) = l, 

At 
Ax 7 

+c(t„,x m _i)] (m = 1, . . . ,M - 1), 
c(i n+ i,x ) = c(i n+1 ,x 1 ), c(i n+ i,x M ) = 1 (n = l,2,...). (20) 



The relative mass gain is given by Eq. (18), where the integral is estimated by means of 
the Euler formula, 



M-1 



$(t n ) = Ax c(t n ,x m ). (21) 

m=0 

The coefficient D is determined by the following procedure. First, we fix some interval 
[Dmin, -Dmax], where the "best-fit" value of D is supposed to be located, and divide this 
interval into I sub-intervals by the points Di = D m i n + iAD (i = 0,1,...,/), where 
AD = (-D max — D min )/I. For any Z)j, Eq. (20) is solved numerically (with the steps 
Ax = 0.05 and At = 0.005 that guarantee the stability of the numerical algorithm), and 
the function $(t n ) is given by Eq. (21). The "best-fit" value of D is found from the 
condition of minimum of the function 



S[ $ «p(**) _ $ num(4) 



2 



tk 

where the sum is calculated over all point, ik, depicted in Figure 1, § cxp (tk) is the relative 
mass gain measured in the test, and <& num (i k ) is determined by Eq. (21). After the "best- 
fit" coefficient, Di, is found, this procedure is repeated for the new interval [-Dj-i, Di+i] 
to ensure an acceptable accuracy of matching observations. 

To show that the numerical results are not affected by our choice of the steps of 
integration, Ax and At, we solve Eqs. (20) and (21) numerically with the "best-fit" 
parameter D and the new steps, Ax = 0.01 and At = 0.0005. No difference is visible 
between the curve plotted in Figure 1 and the curve obtained by numerical integration of 
Eqs. (20) and (21) with the decreased steps in time and the spatial coordinate. 

The above procedure is repeated to match observations in the tests on 4 specimens 
with various thicknesses, and the average diffusivity (over 4 specimens) and its standard 
deviation are determined. 

The same numerical algorithm is employed to fit experimental data on nanocomposites 
with various contents of clay. The difference in the treatment of observations consists in 
the following: 

1. Eq. (20) is replaced by an appropriate finite-difference approximation of Eqs. (9) 
and (11), 
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2. the "best-fit" values are simultaneously determined for 3 adjustable parameters, D, 
C and K. 

Figures 1 to 5 demonstrate fair agreement between the experimental data and the results 
of numerical analysis. 

It is worth noting that our algorithm for the approximation of observations differs 



from that employed in || [11], |29[. In the previous works, the coefficient of diffusion was 
determined by matching the initial parts of the mass gain diagrams (where the function 
$(£) is approximately linear), which did not allow immobilization of water molecules on 
the surfaces of clay layers to be accounted for. 

The average values of D [given D, the diffusivity, D, is determined by Eq. (10)], C 
and K (over 4 tests) are plotted versus the clay content in Figures 11 to 13 (the vertical 
bars stand for the standard deviations). The dependence of diffusivity, D, on the clay 
content, u, is approximated by the phenomeno logical equation 

D = D + D 1 exp(-—\ (22) 

where D , D\ and uq are adjustable parameters. The quantity z/ is determined by the 
steep est- descent method to minimize deviations between the experimental data and their 
approximation by Eq. (22). The coefficients Do and D\ are found by the least-squares 
algorithm. 

The effect of clay content, u, on the steady concentration, C, of immobilized water 
molecules is predicted by the linear function 

C = O, (23) 

where the coefficient C\ is determined by the least-squares technique. 

To describe the influence of clay content, u, on the dimensionless rate, K, of water 
adsorption on the clay platelets, we employ the phenomenological relation similar to Eq. 
(22), 

\ogK = K + K 1 ex P (--), (24) 

where log = log 10 . We use the same value, z/ = 0.55, that was found by fitting observa- 
tions for the coefficient of diffusion with the help of Eq. (22). The material parameters, 
K and Ki, in Eq. (24) are determined by the least-squares method. 

It follows from Eq. (17) that the maximal moisture uptake by the polymeric matrix, 
-Rpoi = Kn°/p, reads 

flpo, = (25) 

Using measurements of the maximal mass gain, R(oo), and the values of C determined by 
matching observations, we calculate _R po i from Eq. (25). This quantity is plotted versus 
clay content, u, in Figure 14. By analogy with Eq. (22), we fit the experimental data by 
the equation 

R pol = R -R ie xp(--), (26) 



where the coefficients Rq and R\ are found by the least-squares algorithm 
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Equations (12) and (17) imply that the maximal moisture uptake by the clay particles, 
-Relay = Knl/p, is given by 

-Relay = -R(oo) — -Rpol- 

We determine -R c i ay from this equality and present this quantity as a function of clay 
content, v, in Figure 15. The experimental data are approximated by the linear function 

-Relay = R\^i (27) 

where the constant R\ is found by the least-squares technique. 

Figures 11 to 15 demonstrate an acceptable agreement between the observations and 
their approximations by phenomenological relations (22) to (24), (26) and (27). 



5 A model for the elastoplastic response 

A nanocomposite is modelled as an equivalent network of chains bridged by junctions. 
Because the viscoelastic behavior of the nanocomposite far below the glass transition 
point is disregarded, the network is treated as permanent (which means that within the 
experimental time-scale, chains cannot separate from junctions as they are agitated by 
thermal fluctuations). The elastoplastic behavior of the network is treated as sliding of 
junctions with respect to their reference positions in the bulk medium. 

For uniaxial deformation, sliding of junctions is determined by a plastic strain e p . 
Adopting the conventional hypothesis that the macro-strain, e, is transmitted to all chains 
in the network by surrounding macromolecules, we find that 

e = e e + ep, (28) 

where e e is an elastic strain. We adopt the mean-field approach, according to which e e 
and e p are treated as average elastic and plastic strains per chain. 

The rate of changes in the strain, e p , with time, t, is assumed to be proportional to 
the rate of changes in the macro-strain e, 

^r(t) = <p(e e (t))j t (t), (29) 

where the coefficient of proportionality, tp, is a function of the elastic strain e e . The 
function p(e c ) vanishes at the zero elastic strain, <p(0) = 0, monotonically increases with 
e e , and approaches some constant a G (0, 1) at relatively large elastic strains. The constant 

a = lim <p(e e ) 

determines the rate of sliding in junctions for a developed plastic flow. 

At small strains, a chain is modelled as a linear elastic solid with the mechanical energy 

1 2 

w = ~pe e , 
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where \i is an average rigidity per chain. Multiplying the energy, w, by the number 
of chains per unit volume, no, we find the strain energy density per unit volume of a 
nanocomposite 



W = -Eel 



(30) 



where E = un is an elastic modulus. 

It is worth noting that E may differ from Young's modulus that is conventionally 
determined as the tangent of the angle between the tangent straightline to a stress- 
strain curve at small strains and the horizontal axis. These two quantities coincide when 
the elastic strain, e e , equals the macro-strain, e. When the stress-strain diagram has a 
pronounced curvature, E exceeds Young's modulus found by the traditional method. 

For isothermal uniaxial deformation, the Clausius-Duhem inequality reads 

„ . . dW . . . . de , . 
Q( t ) =-_(,) + , (i )-( t )>0, 

where Q is internal dissipation per unit volume. Substition of Eqs. (28) to (30) into this 
equality results in 

Q(t) = [a(t) - Ee c (t)(l - V (e c (t)))] ^(f) > 0. 

Assuming the expression in square brackets to vanish, we arrive at the stress-strain rela- 
tion 

a(t) =Ee c (t)[l-<p(e c (t)j\. (31) 
To approximate experimental data, we use the phenomenological equation 



<p(e c ) = a 1 - exp(-— ) 



(32) 



which is determined by two adjustable parameters, a and e > 0. 

Constitutive equations (28), (29), (31) and (32) contain 3 adjustable parameters: 

1. the elastic modulus E, 

2. the rate of developed plastic flow a, 

3. the strain, e, that characterizes transition to a steady plastic flow. 

These quantities are found by matching the experimental data depicted in Figures 6 to 
10. 



6 Fitting observations in mechanical tests 

It follows from Eqs. (28), (29), (31) and (32) that in a uniaxial tensile test, the longitudinal 
stress, a, is given by 



a(e) = E(e - e p ){l - a[l - exp(-^-^)] }, 



(33) 
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where the plastic strain, e p , satisfies the nonlinear differential equation 



de 



1-exp p -) , e p (0)=0. (34) 



We begin with matching the experimental data in a test on the neat vinyl ester resin. 
To find the constants, E, a and e, we fix some intervals [0,a max ] and [0,e max ], where the 
"best-fit" parameters a and e are assumed to be located, and divide these intervals into 
/ subintervals by the points a, = iAa and Ej = jAe = 1, . . . , I) with Aa = a max /J 
and Ae = e mBX /I. For any pair, {cti,£j}, Eq. (34) is integrated numerically by the 
Runge-Kutta method with the step Ae = 1.0 • 10~ 5 . Given {di,Ej}, the elastic modulus 
E = E(i,j) is found by the least-squares technique from the condition of minimum of the 
function 2 

F = ^2[a cxp (e k ) - a num (e k ) , 

where the sum is calculated over all experimental points, e k , depicted in Figure 6, cr cxp 
is the stress measured in a tensile test, and a num is given by Eq. (33). The "best-fit" 
parameters a and e minimize the function F on the set {a,, Ej — I)}. Fitting 

the observations results in e = 2.11 • 1CT 3 . 

We fix this value of e and proceed with matching observations in tensile tests on other 
specimens by using only two material constants, E and a, which are determined by the 
above algorithm. Figures 6 to 10 demonstrate fair agreement between the observations 
and the results of numerical simulation. 

The average values of E and a (over 4 samples for any concentration of clay) are 
depicted in Figures 16 and 17, where the vertical bars stand for the standard deviations. 
Figure 16 shows that E is practically independent of the clay content. 

The effect of clay content, z/, on the rate of developed plastic flow, a, is described by 
the phenomenological equation 

a = a + axexpf V (35) 

v *V 

where vq is found by matching experimental data for the coefficient of diffusion, D, and the 
coefficions, a and ai are determined by the least-squares method. Figure 17 demonstrates 
good agreement between the observations and their approximation by Eq. (35). 



7 Discussion 

Changes in the elastic modulus of a hybrid nanocomposite may be explained by the influ- 
ence of two factors. On the one hand, it is conventionally accepted that reinforcement of 
a polymeric matrix by rigid particles results in an increase in the elastic modulus of the 
composite. On the other hand, clay platelets in a melt-intercalated nanocomposite screen 
the macro-strain, which is not transmitted to chains in the close vicinity of nanoparticles 
(in particular, to chains intercalated into the galleries between silica sheets). As a result, 
meso-regions of "occluded" polymer arise similar to the domains of occluded rubber in 
particle- reinforced elastomers |3l| . These occluded regions are not deformed, which im- 



plies that the total mechanical energy of a nanocomposite (per unit volume) and, as a 
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consequence, its elastic modulus, E, decrease with v. Our observation that the elastic 
modulus, E, is independent of the clay content, is, (Figure 16) may be ascribed to a 
cumulative effect of these two micro-mechanisms. 

Weak changes in the elastic modulus of intercalated nanocomposites and even its de- 
crease with clay content were observed by several authors. Bharadwaj et al. |§ and Xu et 



al. P2| reported a decrease in the elastic modulus of nanocomposites with polyester and 
epoxy resin matrices, respectively, and associated this behavior with a strong decay in the 
degree of cross-linking induced by the presence of clay particles [§j . Lepoittevin et al. 
observed a decrease in Young's modulus for poly(e-caprolactone) filled with a non-modified 
MMT clay and a propounced growth of the elastic modulus of the nanocomposites with 
organo-modified clay particles, but did not provide explanations for their findings. Kim et 
al. suggested that changes in elastic moduli and toughness of nanocomposites may be 
attributed to nucleation of micro-voids in the vicinity of stacks of silicate sheets that are 
thought of as initiation sites for cavitation. This indicates that physical mechanisms gov- 
erning the effect of clay concentration on elastic moduli of hyrid nanocomposites remain 
a subject of debate. 

Figure 17 demonstrates that the rate of developed plastic flow, a, monotonically de- 
creases with the clay content, v. This trend is in agreement with a conventional standpoint 
[ 35| that the presence of nanoparticles in an intercalated nanocomposite leads to reduction 
of segmental mobility of polymeric chains confined to galleries between clay platelets . It 
is also confirmed by DSC measurements |29| that show a pronounced increase in the glass 
transition temperature, T g , with the clay content (the growth of T g confirms a decrease 
in molecular mobility of the host matrix). 

Figure 11 reveals a similar decrease in the rate of transport of moisture molecules 
through the nanocomposite. Fair approximation of the experimental data depicted in 
Figures 11 and 17 by the exponential functions with a fixed constant, z/ , implies that the 
same micro-mechanism may be responsible for the effect of clay content on the quantities 
D and a. According to the free volume concept |B6j, B7J, an elementary mass-transport 



event is rationalized as follows. At some instant, t, a small penetrant molecule in a host 
polymer is at rest (thermal oscillations are disregarded) waiting for creation of a micro-hole 
nearby. When a hole with a sufficient size is formed (due to thermal motion of chains), the 
penetrant molecule hops into the hole, and remains in the new position until a new hole 
arises in its neighborhood. Although this picture is oversimplified (such important factors 
are neglected as clustering of water molecules, their weak bonding to polymeric chains, 
obstacles on the path of a penetrant molecule caused by the presence of inclusions), 
it provides a clear indication that mass transport is strongly influenced by mobility of 
chains. Resemblance of the graphs presented in Figures 11 and 17 implies that moisture 
diffusion and elastoplasticity of a nanocomposite may reflect the same phenomena at the 
micro-level associated with molecular mobility of the polymeric matrix. 

According to the above scenario of the mass-transport process as a sequence of hopping 
and waiting events, slowing down of water diffusion (observed as a pronounced decrease 
in diffusivity D) implies that the waiting times strongly increase. As a consequence, 
the probability of aggregation of penetrant molecules grows, which results in formation 
of small clusters of water in a nanocomposite. Clustering of small molecules in solid 



polymers was observed in several studies, see |3q, 39, 40] and the references therein. 
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According to the free volume theory, clustering of water molecules implies an increase 
in the total mass uptake by the polymeric matrix (because penetrant molecules weakly 
bounded into aggregates occupy less volume than individual molecules). This conclusion 
is confirmed by the experimental data for R po \ plotted in Figure 14. It follows from 
phenomenological equations (22) and (26) that an increase in the moisture uptake with 
clay content strongly correlates with a decrease in diffusivity (the same parameter, u , is 
used to fit the observations depicted in Figures 11 and 14). 

The coefficient of diffusion, D, may serve as an average (over a macro- volume) rate of 
molecular transport. The local rate of moisture diffusion can substantially differ in the 
close vicinities of the clay particles (where molecular mobility is severely reduced due to 
intercalation of chains into galleries between platelets) and in the bulk (where mobility of 
chains is weakly affected by the presence of filler). This implies that the mass transport 
in a nanocomposite becomes substantially heterogeneous. 

According to the concept of water clustering, noticeable slowing down of diffusion near 
inclusions implies that clusters of penetrant molecules are mainly formed in neighborhoods 
of clay particles. This means that immobilization of a water molecule on the hydrophilic 
surface of a clay layer may occur in two different ways: 

1. as an one-step event, when an individual molecule merges with an unoccupied site 
of the surface of a silicate sheet, 

2. as a two-step event, when a water molecule (i) separates from a cluster, and (ii) 
merges with an appropriate site. 

At small concentrations of clay, when diffusivity, D, is rather large and clustering of water 
molecules is negligible, immobilization of small molecules takes place by the one-step 
mechanism, which means that the rate of bounding of penetrant molecules is relatively 
high. With an increase in the clay content (which is tantamount to a decrease in the 
average diffusivity, D, see Figure 11), the one- step mechanism of bounding is replaced 
by the two-step mechanism. According to the proposed scenario, these changes in the 
immobilization mechanism result in a noticeable deceleration of the process of water 
bounding to silica sheets, because separation of a water molecule from a cluster requires 
an additional energy provided by thermal fluctuations. This slowing down (reflected 
by the model as a decrease in the coefficient K) is confirmed by the experimental data 
presented in Figure 13. This figure demonstrates that K is reduced by at least two orders 
of magnitude with the growth of clay content. 

It is natural to assume that the total moisture uptake by the clay particles, R c \ ay , is 
proportional to their concentration, v. The parameter -R c ia y is plotted in Figure 15 versus 
clay content. The experimental data show that the linear dependence (27) is in good 
agreement (within the experimental uncertainties) with the observations. This conclusion 
may be treated as a confirmation of our kinetic equations for anomalous moisture diffusion. 

Another confirmation of the model is provided by the fact that the diffusivity, D, 
becomes practically independent of the clay content when v exceeds 5 wt.-% (Figure 
11). This conclusion is in accord with numerous observations that reveal substantial 
improvement of material properties of hybrid nanocomposites with clay content in the 
range of several (conventionally, 5) wt.-%, whereas at higher clay fractions, material 
parameters are weakly affected by the MMT concentration. 
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8 Concluding remarks 



Diffusion tests and uniaxial tensile tests with a constant strain rate have been performed 
on an intercalated nanocomposite with vinyl ester resin matrix and montmorillonite clay 
filler at room temperature. Observations in diffusion tests show that moisture diffusion 
in the neat resin is Fickian, whereas it becomes noticeably anomalous (non-Fickian) with 
the growth of the clay content. This transition is attributed to immobilization of the pen- 
etrant molecules on hydrophilic surfaces of clay layers. Experimental data in mechanical 
tests demonstrate that the response of the neat vinyl ester resin is strongly elastoplastic, 
whereas plastic strains strongly decrease with the clay content. This observations are 
ascribed to slowing down of molecular mobility in the host matrix driven by confinement 
of chains in galleries between clay layers. 

Constitutive equations have been developed for anomalous moisture diffusion through 
and for the elastoplastic behavior of a nanocomposite. These relations are determined by 3 
adjustable parameters (for each process) that are found by fitting the experimental data. 
Fair agreement is demonstrated between the observations and the results of numerical 
simulation. 

The following conclusions are drawn: 

1. The diffusivity, D, and the rate of developed plastic flow, a, decrease with clay 
content following similar dependences. The effect of clay concentration, v, on these 
quantities may be explained by slowing down of molecular mobility in the polymeric 
matrix. 

2. The moisture uptake by the host matrix grows with clay content. With reference 
to the free volume concept, this observation is explained by clustering of water 
molecules in the close vicinity of stacks of platelets (where diffusivity falls down 
dramatically). 

3. The rate of immobilization of water molecules on the hydrophilic surfaces of inor- 
ganic sheets decreases noticeably with v, which is associated with transition from 
the one-step to the two-step mechanism of immobilization. 

An important advantage of our model compared to the tortuous path concept is that 
it can explain (in a unified manner) a decay in diffusivity of a hybrid nanocomposite 
(with the clay content) together with an increase in the maximal moisture uptake by the 
polymeric matrix and a reduction of its elastoplastic properties. 

It is worth noting that the correlations revealed between a decrease in diffusivity and 
a decrease in the rate of a steady plastic flow provide a way to substantially reduce 
the number of (time-consuming) diffusion tests by replacing a part of experiments by 
mechanical tests that do not require long-time observations. 
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List of figures 



Figure 1: The relative water uptake $ versus the reduced time i (h^/mm). Circles: 
experimental data (21 = 0.49953 mm, v = 0.0 wt.-%). Solid line: results of numerical 
simulation 

Figure 2: The relative water uptake <£> versus the reduced time t (ha /mm). Circles: 
experimental data (21 = 0.17318 mm, v = 0.5 wt.-%). Solid line: results of numerical 
simulation 

Figure 3: The relative water uptake $ versus the reduced time t (ha" /mm). Circles: 
experimental data (21 = 0.18344 mm, v = 1.0 wt.-%). Solid line: results of numerical 
simulation 

Figure 4: The relative water uptake $ versus the reduced time t (li^/mm). Circles: 
experimental data (21 = 0.1397 mm, v = 2.5 wt.-%). Solid line: results of numerical 
simulation 

Figure 5: The relative water uptake $ versus the reduced time i (ha /mm). Circles: 
experimental data (21 = 0.18452 mm, v = 5.0 wt.-%). Solid line: results of numerical 
simulation 

Figure 6: The stress a MPa versus strain e. Circles: experimental data (v = 0.0 wt.-%). 
Solid line: results of numerical simulation 

Figure 7: The stress o MPa versus strain e. Circles: experimental data (v = 0.5 wt.-%). 
Solid line: results of numerical simulation 

Figure 8: The stress a MPa versus strain e. Circles: experimental data (v = 1.0 wt.-%). 
Solid line: results of numerical simulation 

Figure 9: The stress a MPa versus strain e. Circles: experimental data (v = 2.5 wt.-%). 
Solid line: results of numerical simulation 

Figure 10: The stress a MPa versus strain e. Circles: experimental data (v = 5.0 wt.-%). 
Solid line: results of numerical simulation 

Figure 11: Diffusivity mm 2 /s versus the clay content v wt.-%. Circles: treatment 

of observations. Solid line: approximation of the experimental data by Eq. (22) with 
D = 3.53 and D 1 = 5.93 

Figure 12: The dimensionless parameter C versus the clay content v wt.-%. Circles: 
treatment of observations. Solid line: approximation of the experimental data by Eq. 
(23) with Ci = 0.057 

Figure 13: The rate of adhesion K versus the clay content v wt.-%. Circles: treatment 
of observations. Solid line: approximation of the experimental data by Eq. (24) with 
K x = 0.29 and K x = 4.80 

Figure 14: The moisture adsorption by the polymer i? pol % versus the clay content v 
wt.-%. Circles: treatment of observations. Solid line: approximation of the experimental 
data by Eq. (26) with R = 0.90 and R x = 0.49 

Figure 15: The maximal moisture uptake by the clay particles -R c ia y = Kn°/p % versus 
the clay content v wt.-%. Circles: treatment of observations. Solid line: approximation 
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of the experimental data by Eq. (27) with R± = 0.053 

Figure 16: The elastic modulus E GPa versus the clay concentration v wt.-%. Circles: 
treatment of observations. Solid line: approximation of the experimental data by the 
constant E = 4.60 GPa 

Figure 17: The rate of developed plastic flow a versus the clay content v wt.-%. Circles: 
treatment of observations. Solid line: approximation of the experimental data by Eq. 
(35) with a = 0.106 and ai = 0.106 



20 




70.0 



70.0 



Figure 2: 



21 




70.0 



70.0 



Figure 4: 



22 




70.0 



Figure 5: 



40.0 



a 



0.0 




0.0 e 0.01 



Figure 6: 



23 




Figure 8: 



24 




Figure 10: 



25 



12.0 





27 




Figure 15: 




0.0 



v 5.0 



Figure 16: 



28 



0.3 




Figure 17: 



29 



